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Abstract. We study a magnitude of possible over /underestimation of 
the actual 7-ray luminosity ^actual of a millisecond pulsar when using so- 
called pseudo luminosity Lpseudo which is inferred from a phase-averaged 
flux. Both, -^actual and ^pseudo are calculated with the numerical code 
that models the millisecond pulsar magnetospheric activity in the space 
charge limited flow approximation with unscreened accelerating electric 
field. The behaviour of -^^pseudo/-^^ actual depending on viewing angle is 
analysed in the energy bands corresponding to the operational energy 
bands of Fermi GST and H.E.S.S. II. 

1. Introduction 

Magnetospheric activity of the millisecond pulsars is numerically modelled in 
order to test the correctness of pulsar theory and also in order to develop the 
tool that would predict the detectability of these pulsars with high energy ob- 
servatories. 

Our main reason for carrying out the numerical modelling of radiative pro- 
cesses taking place in millisecond pulsar magnetosphere is to check the accuracy 
of the so-called pseudo luminosity Lpseudo which is inferred from the observed 
phase-averaged flux as an estimation of the pulsar luminosity in 7-rays. In the 
calculations the unscreened electric field that includes general relativity effects 
is taken into account. The calculations are made for a typical millisecond pulsar 
with the spin period P = 2.3 ms and the magnetic field strength at the pole 
Bpc = TG. We restrict our analysis to the energy range of 7-rays where 

the emission is due to curvature process. Analogous treatment of X-ray char- 
acteristics of millisecond pulsars would require taking into account synchrotron 
radiation due to secondary pairs as well as thermal component and this is beyond 
the scope of this work. 

The main features of the numerical model are presented in Section 2; Sect. 
3 shows the preliminary results of the modelling. The concluding remarks are 
gathered in Sect. 4. 



2. Numerical Model 

The radiative processes taking place in pulsar magnetosphere are modelled in 
3D. It is assumed that charged particles move along open magnetic field lines 
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Figure 1. The accelerating electric field for a millisecond pulsar charac- 
terised with P = 2.3 ms, Spc = lO'^ TG, a = 20° and ^ = 0.7. The total 
electric field | (solid line) is a combination of a formula describing the elec- 
tric field near the neutron star i?noar (dashed line) up to a distance h ~ 0.7 
of the polar cap radius where a formula depicted with i?far (dot-dashed line) 
becomes applicable. 



where their flow is tr eated within the framework of the space charge limited flow 
(jFawley et al.l . Il977l ). The particles accelerate in the unscreened electric field 
(Muslimov Sz Harding, 1997; Harding & Muslimov, 1998) tha t includes the gen- 
eral r elativistic efl'ect of dragging of inertial frames of reference ( Muslimov fc TsyganI . 
I1992I ). It is predicted that for the majority of millisecond pulsars the acceleration 
voltage drop that develop s above polar cap is not limited by the development 
of a pair formation front ( Harding et al.l . l2002l ). Thus, particles are accelerated 
even at high altitudes and it is justifled to describe the accelerating electric field 
as being unscreened. 

In the magnetospheres of millisecond pulsars a variety of radiative processes 
may take place. In this work, however, we restrict the considered processes to 
curvature radiation and the magnetic absorption of curvature photons without 
tracing the created secondary particles because it is the curvature emission which 
(according to polar ca p models) domina tes completely the 7-ray energy range 
in millisecond pulsars (|Bulik et al.l . I200Q| ). The escape rate of primary electrons 
from the neutron star surf ace is assumed to be the Go ldreich-Julian rate in the 
general relativistic regime ( Harding fc Muslimovl . [19981 ) . 

Figure H] demonstrates the behaviour of the electric field - the component 
parallel to the local magnetic field - above the polar cap (solid line). The accel- 
erating field is a combination of the formula describing the electric field in the 
vicinity of a pulsar (dashed line; Dyks & Rudak, 2000) and the formula applica- 
ble to distances beyond ~ 0.7 of the polar cap radius (dot-dashed line; Muslimov 
& Harding, 1997). In Fig.Hlthe case of pulsar with P = 2.3 ms, Bpc = 10"^ TG, 
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a = 20° and = 0.7 is shown. Here a is an inchnation angle of the magnetic 
axis with respect to the spin axis and ^ = where is a magnetic colatitude 

and 9{'q) is a half-opening angle of the polar magnetic flux tube; r/ = 1 + -^^j 
where h is an altitude above the neutron star surface. 



3. Results 

The parameters of the modelled typical millisecond pulsar that were used in 
calculations are: radius i^NS = 10^ cm, mass Mns = 1-4 Mq , moment of 
inertia Ins = 10^^ g cm^, spin period P = 2.3 ms, magnetic field strength at 
the pole Bpc = 10~^ TG, inclination angle a = 20° and spin down luminosity 
Lsd ^ ^TT^c-^B^^ i2^s P-^ ~ 5.2 • 1035 erg s'^ 

We present the results for a single millisecond pulsar in the form of pho- 
ton density maps, which demonstrate the directional dependence of the photon 
emission in a chosen energy band. Such map is a contour plot where the photon 
density defined as 

Je^.^ dt dE dC dip 

is coded in grey scale (A^ is a number of photons). It is presented in the coordinate 
system where the pulsar rotation phase (j) = ip/2'K \s on the x-axis and the 
observing angle C is on the y-axis. Note that these m aps are present e d in a 
different manner than the p ho ton densi t y ma ps in e.g. iDvks Rud"^ (120031 1. 
Frackowiak fc RudakI (l2005l ) or lHardind (l2007l ). 



Two energy bands have been chosen in which the photon density maps of the 
modelled millisecond pulsar are presented (Fig. [2]). The low energy band stretches 
from 100 MeV to 30 GeV, while the high energy band covers the energies from 30 
GeV up to 10 TeV. The energy division corresponds to the energy bands of Fermi 
GST (the low energy band) and H.E.S.S. II (the high energy band). Photons 
with energies below 30 GeV (Fig. O top panel) are emitted in a range of the 
viewing angle C, that is wider in comparison with the range for the high energy 
band (Fig. [21 bottom panel). The emission cone of the high energy photons is 
centred at the pulsar magnetic axis. This centring of emission is also traceable in 
the horizontal cuts made through the photon density maps for chosen observers 
(Fig. [31). In the energies above 30 GeV the level of the photon density drops 
dramatically if the observer moves away from the magnetic axis (Fig. [3l right 
panel). The decrease in the photon density level is also present in the lower 
energy band (Fig. [3l left panel); however, it is not so rapid. This behaviour is 
clearly visible if one compares the light curves for C = 5°, 20° and 35° in different 
energies. The overall photon density level is lower for the high energy photons. 

To have an idea how accurate the estimations of the average power radiated 
by the pulsar are, we have calculated a pseudo luminosity Lpseudo for both energy 
ranges. -Lpseudo emulates the power inferred by an observer viewing the pulsar 
at an angle Q. The emitted radiation is confined to a solid angle f^pseudo = 1 sr. 
The formula for the pseudo luminosity in the energy range from £^min to E'n 
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Figure 2. Photon density maps -/V<^,i^(£^min, ^^max) in the energy ranges: 
from 100 MeV to 30 GeV (top panel) and from 30 GeV to 10 TeV (bottom 
panel). The pulsar rotation phase (j) = if/2TT is on the x-axis and the 
viewing angle C in degrees is on the y-axis. White horizontal lines indicate 
observers for which corresponding light curves have been plotted in Fig. [3l 
( = 5°, 20°, 35°, 50° are depicted with soHd, dotted, dashed and dot-dashed 
Hne, respectively. See discussion of the photon density maps in the text. 
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Figure 3. Pulsar light curves for the observers with C = 5°, 20°, 35°, 50° 
in the energy range from 100 MeV to 30GeV (left panel) and from 30 GeV 
up to 10 TeV (right panel). The light curves are due to the horizontal cuts 

through the photon density maps A^c,i^(i?min, -E-max) presented in Fig. [21 The 
values of the photon density are in arbitrary units. 



for a viewing angle ( is thus given as 

^pscudo Z"'^'"'"' /■^'^ dN 

-^pseudo(C; -E'min, -E'max) = „ ° / / .a . dip dE . (2) 

27r smC Jo dt dE dC dip 

-^^pseudo is a theoretical counterpart of Lobs inferred from observations in the given 
energy range {Emin,Emax), which can be described with the following formula: 

-^obs (-^mim -^'max) — ^arb ' ^ ' (-^(-E'min) -^max)) ) (3) 

where Oarb is a solid angle chosen arbitrarily (1 steradian or sometimes 2tt stera- 
dians), D is a distance to the object and (-F(Sjnin, -Emax)) is a phase-averaged 
energy flux within (-Emin, -Emax) inferred from observations. This formula is fre- 
quently used to asses the pulsar luminosity from the observed phase-averaged 
flux and the known distance to the object. The modelled pseudo luminosities 
-^pseudo will therefore test the accuracy of the observed luminosities Lobs versus 
the actual luminosity Lactual- 

The theoretical pseudo luminosity spectra are presented in the top panel 
of Fig. m In the same figure the spectrum of the actual power radiated by the 
pulsar from one pole is shown with solid line. 

For each observer we can construct the bias factor 

^_A-eudo ^ 
Lactual 

where Lpseudo is given by Eq. [2] and Lactual is defined as 

^actual (ii^min, i?„.ax) = E ^^^^^ dip dE . (5) 
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Figure 4. Top panel: Pseiido luminosity spectra iE,pscudo:C the ob- 
servers with ( = 5°, 20°, 35°, 50° are shown. Actual output from each pole 
^E, actual is presented here with solid line. Here Le = j^; iE,pseudo;C = 

fet E and LE.ac.uai = /.^ I^f ^ ^of^ '^C d^. 

Bottom panel: The dependence of the bias factor b on the viewing angle ( 
scaled with the inclination angle a is presented in the energy range from 100 
MeV to 30 GeV and from 30 GeV to 10 TcV with a solid and a dot-dashed 
line, respectively. See discussion of the bias factor in the text. 
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Both i^pseudo and -^^actuai are calculated within the same energy range (-Emin, 
-E'max)- The dependence of the bias factor on the viewing angle is presented in the 
bottom panel of Fig. [H In the low energy range (solid line) the bias factor equals 
1 for C — 0.8 a. The observers with ( < 0.8 a tend to overestimate the pulsar 
actual luminosity (Eq. [5|) using the pseudo luminosity (Eq. [2]) , while observers 
with C ^ 0.8 a underestimate -^^actuai- In the high energy range (dot-dashed line) 
the situation is different due to the opening angle of radiation being smaller 
than the inclination angle a and b = 1 for ^ ~ 0.25 a and C — 1.2 a. Viewing 
the pulsar with C in the range from 0.25 a to 1.2 a results in overestimation of 
Lactual- However, the overestimation is rather insignificant {b is less than 1.5 in 
maximum). On the other hand, viewing it with C lower than 0.25a or higher 
than 1.2 a leads to underestimation of the pulsar actual luminosity. 
Finally, we calculate the gamma-ray efficiency factoiQ : 

^7 = —P • (6) 

Here, -^actual (Eq. [5]) is calculated within the energy range stretching from 100 
MeV up to 10 TeV. This factor demonstrates what part of the pulsar spin-down 
luminosity Lgd is converted into the gamma radiation. The gamma-ray efficiency 
factor for the modelled millisec ond pulsar is r? ^ ^ 0. 05. This result is in good 
agreement with rj^ estimated by iHarding et al.l ( 20021 ). 



4. Conclusions 

The preliminary results of the numerical modelling of millisecond pulsars have 
been presented in this paper. We have calculated the photon density maps in 
two energy regimes corresponding to the operational energy bands of Fermi GST 
and H.E.S.S. H. We have also computed the pseudo luminosities for different 
viewing angles and on their basis the plot presenting the behaviour of the bias 
factor has been constructed. From the comparison of the bias factor - viewing 
angle dependence in different energy bands we may infer information on accuracy 
of the so-called pseudo luminosity -Lpseudo as an estimation of the pulsar actual 
luminosity ^actual in 7-rays. It turns out that there are very few favourable (" for 
which -Lpseudo is a good estimate of -^actual- Most often Lactual tends to be un- 
der/overestimated. The gamma-ray efficiency factor for the modelled millisecond 
p ulsar has been c alcula ted. Its value r/^ ~ 0.05 is in agreement with the values 
in lHarding et al. I (l2002l ). However, further work concerning modelling, especially 



the inclusion of radiative processes like synchrotron radiation and inverse Comp- 
ton scattering important for pulsar X-ray emission, has to be done in order to 
draw accurate conclusions about millisecond pulsars. Also simulations for other 
sets of pulsar parameters, especially for different inclination angles a, have to be 
carried out. 

We intend to use the results of modelling as a tool for predicting the de- 
tectability of millisecond pulsars with high energy observatories like Fermi GST, 



^i'actuai (Eq. [U is the power output only from one pole. Thus, in order to compare it with Lsd 
we have to multiply it by 2 to include the output from both poles. 
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H.E.S.S. II and even CTA (iFrackowiak Rndakl . EiORl V Furthermore, we plan 
to develop the population synthesis studies of high energy pulsars that would 
incorporate the results of the modelling described briefly in this paper. 
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